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Report  to  Director  of  Research  Materials  2  under  former  Ministry  of  Technology 

Research  Agreement  No.PD/57/021 


UNLIMITED 


FORWORD 


This  report  was  prepared  by  the  Department  of  Engineering  Science , 
Oxford  University  for  the  Director  of  Research  Materials  2, 
Procurement  Executive,  Ministry  of  Defence  under  former  Ministry 
of  Technology  Research  Agreement  No. PD/57/021.  The  work  was 
conmenced  in  October,  1966.  The  work  was  planned,  progressed  and 
administered  by  Mat(0)l  under  the  direction  of  ADR/Mat2(0) , 
I-rocurement  Executive,  Ministry  of  Defence.  The  work  was  performed 
by  Mr.  R.  G.  C.  Arridge  assisted  by  Mr.  R.  C.  Stone  under  the 
supervision  of  Dr.  N.  G.  McCrum. 

This  report  covers  work  performed  between  1966  and  1969. 
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SUMMARY 

V  ~ 

Thve  report  describes  the  measurement  of  shear  modulus  and  damping  of 
epoxy  resin  samples  undirect ionally  reinforced  with  silica  fibres.  The 
measurements  have- been  performed  over  a  temperature  range  including  the 
glass-rubber  transition  of  the  resin  on  samples  cut  with  fibres  at  various 
angles  to  the  specimen  axis.  Fibre  concentrations  up  to  74Z  have  been  used. 

The  results  have  been  compared  with  theoretical  values  derived  from 
three  different  theories  and  fair  agreement  found  for  the  glassy  state  of 
the  resin.  Anomalies  wer?  found  in  the  rubber  modulus  and  in  the  shape  and 
position  of  the  main  absorption  peak  for  the  resin.  The  causes  of  these 
anomalies  are  discussed  and  satisfactory  explanations  put  forward. 

Section  3  of  the  report  presents  a  theory  of  fracture  of  unidirectional 
fibre  reinforced  materials  in  which  the  rate  at  which  elastic  energy, 
released  by  fibre  breakage,  can  be  absorbed,  is  calculated. 

Analogies  of  fibre  reinforcement  theory  and  chain  fracture  theories  in 
polymers  are  discussed. 
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1.  INTRODUCTION 


The  improvement  of  the  properties  of  synthetic  resins,  of  both  the  thermosetting 
and  the  thermoplastic  types,  by  the  incorporation  of  high  strength,  high 
modulus  fibres  is  by  now  a  well-known  technique.  A  great  deal  of  practical 
knowledge  on  glass  reinforced  polyester  and  epoxide  resins  has  been  obtained  and 
some  theoretical  studies  of  the  elastic  moduli  and  the  fracture  mechanics  of 
such  assemblies  have  also  been  made,  (see  e.g.  Broutman-Krock  (20)), 

The  number  of  variations  possible  in  a  fibre  reinforced  material  allows  the 
engineer  to  consider  designs  in  which  the  most  efficient  use  of  material  can 
be  planned  for  any  particular  component.  However  before  any  exact  theory  can 
be  established  the  properties  of  the  "unit  cell'1  -  the  fibre-matrix  unit  must 
be  determined  and  understood  in  terms  of  the  mechanical  properties  of  each 
component . 

In  the  general  case  neither  the  fibre  nor  t.he  matrix  material  is  an  elastic 
solid  in  the  Hookean  sense.  The  fibre  usually  approximate;:  better  to  this 
ideal  than  does  the  matrix  although  carbon  fibres  show  a  marked  anisotropy  of 
elastic  moduli  and  some  other  types  of  fibre  show  viscoelastic  properties.  A 
resin  matrix  is,  to  a  first  approximation,  linearly  viscoelastic  while  a  metal 
one  possesses  a  yield  point  and  the  characteristics  of  a  crystalline  material 
with  defects. 

Several  theories  have  now  been  put  forward  for  predicting  the  elastic  behaviour 
of  a  fibre  composite  from  its  constituents.  Hashin  and  Rosen,  U)  and  Hill,  (2) 
give  upper  and  lower  bounds  on  the  elastic  moduli  which  are  rigorous  for  the 
models  used.  Whitney  and  Riley  (3)  solve  the  problem  of  the  composite  cylinder 
and  use  this  model  to  predict  properties  of  an  assemblage  of  such  cylinders. 

Savin  and  Van  Fo  Fy  W  solve  the  equations  of  equilibrium  for  a  hexagonal 
array  of  fibres  in  a  matrix  under  applied  stress  using  an  elliptic  function 
expansion  and  integrate  to  find  the  composite  elastic  moduli.  A  numerical 
solution  of  the  stress  equations  for  a  square  array  has  also  been  produced  by 
Heaton  (5)  and  integrated  to  yield  the  elastic  moduli.  Similar  calculations 
have  also  been  made  by  Pickett  (6)  and  Chen  and  Cheng. A  very  full  review 
of  the  theories  is  given  by  Chamis  and  Sendeckyj . (8) 

There  has  been  little  experimental  confirmation  of  any  of  the  theories  except 
for  the  case  of  longitudinally  applied  stresses.  The  results  of  Whitney  and 
Riley  for  example  for  shear  stress  indicate  that  a  large  discrepancy  may  exist 
between  predicted  values  and  those  actually  measured. 

In  a  specimen  reinforced  with  fibres  all  parallel  to  its  longitudinal  axis  the 
longitudinal  Young's  modulus  is  greatly  affected  by  reinforcement.  Seventy 
percent  glass  reinforcement  increases  the  modulus  14  times  and  a  similar 
proportion  of  carbon  or  boron  fibres  increases  it  70  times,  Yet  the  transverse 
Young's  modulus  and  the  longitudinal  shear  modulus  are  only  increased  at  the 
most  by  a  factor  of  4  or  5  over  those  of  the  pure  resin.  Consequently  the 
composite  elastic  properties  most  affected  by  the  properties  of  the  resin  are 
the  transverse  Young's  modulus  and  the  shear  moduli. 

There  appear  to  be  no  published  results  on  fibre  reinforced  viscoelastic 
materials  which  can  be  related  to  theory  when  the  time  dependent  quantities 
are  substituted  in  the  theoretical  equations. 

The  work  reported  was  therefore  undertaken  in  an  attempt  to  discover  whether 
the  relations  derived  by  Hashin,  Hill,  Savin  and  Van  Fo  Fy  and  others  were 
applicable  to  fibre  reinforced  systems  with  time  dependent  mechanical  properties. 


2.  PROGRAMME 

In  order  to  obtain  a  high  modulus  in  the  fibre  direction,  commercially  produced 
fibre-resin  composites  incorporate  as  many  fibres  as  possible.  Consequently, 
fibres  are  packed  in  bundles  and  the  array  cured  under  pressure  in  order  to 
eliminate  voids  and  expel  excess  resin.  The  resulting  material  usually  has  an 
irregular  distribution  of  fibres  and  in  the  fibre  direction,  numerous  cross¬ 
overs  must  occur.  Although  auch  a  material  can  be  considered  as  approximating 
to  a  hexagonal  close  packed  array,  it  is  likely  that  important  physical 
differences  may  exist  because  of  the  irregularity  of  the  packing.  Furthermore, 
the  properties  of  the  fibres  and  resin  used  in  the  composites  are  not  disclosed. 

The  object  of  this  research  was  to  make  a  study  of  fibre  concentration  and  its 
effect  on  elastic  properties.  It  was  therefore  necessary  to  experiment  with 
composites  with  accurately  spaced  fibres  in  a  regular  array.  The  use  of  thick 
rods  which  could  easily  be  arranged  in  any  desired  array,  was  ruled  out  on  the 
grounds  that  to  achieve  a  macroscopic  homogenity,  very  large  specimens  would  be 
required.  This  would  mean  that  the  forces  needed  to  produce  measurable  deflections 
and  thus  enable  moduli  to  be  calculated,  would  be  large  and  difficult  to  obtain 
in  the  laboratory.  In  addition,  there  would  be  doubt  as  to  the  applicability  of 
the  results  to  specimens  with  fibres  of  two  or  more  orders  of  magnitude  smaller. 

It  was  therefore  decided  to  develop  a  method  of  making  ordered  arrays  of  fibres 
of  any  desired  concentration  which  could  be  infiltrated  with  resin  without  the 
use  of  pressure. 

More  than  45  laminates  with  unidirectional  reinforcement  were  prepared  and 
examined.  Specimens  were  machined  from  each  laminate  with  axes  at  different 
angles  to  the  fibres  and  their  torsional  moduli  and  damping  factor  determined 
using  a  torsional  pendulum.  The  determinations  were  made  over  a  temperature 
range  which  included  the  glass  transition  point  of  the  resin  (65°C).  After  the 
mechanical  properties  had  been  measured,  the  specimens  were  cut  up;  fibre 
concentrations  were  determined  on  the  cylindrical  centre  section  and  a  microscopic 
examination  made  on  the  residue.  The  results  obtained  were  then  compared  with 
those  predicted  by  using  the  formulae  postulated  by  three  previous  workers  in  the 
field. 


3.  WORK  CARRIED  OUT  AND  RESULTS  OBTAINED 

3.1.  Manufacture  of  the  fibres  for  reinforcement 

Fibres  of  pure  vitreous  silica  approximately  90p  in  diameter  were  drawn 
from  3  mm  diameter  rods  by  heating  in  an  oxygen/coal  as  flame.  If  a 
finish  was  required  on  the  fibre  surface,  it  was  passed  through  a 
solution  of  10Z  A1100  in  distilled  water.  (A1100  is  the  Union  Carbide 
code  for  y-aminopropyl  triethoxy  silane.) 

3.2.  Preparing  the  fibre  mat 

The  fibres  were  wound  on  to  an  accurately  turned  cylinder  on  a  lathe  using 
the  lead  screw  to  space  the  filaments  by  the  required  amount.  The  cylinder 
was  150  mm  in  diameter  and  300  mm  long,  enabling  four  bands  of  fibres  of 
about  65  mm  width  to  be  wound  on  it.  Longitudinal  strips  of  25  mm  wide 
adhesive  tape  were  secured,  sticky  side  up,  to  the  drum  in  three  places 
around  its  circumference  before  winding  so  ‘■.hat  the  fibres  could  stick  to 
it  and  retain  their  lateral  displacement.  After  completing  the  winding  a 
second  strip  of  adhesive  tape  sticky  side  down  was  placed  over  the  first 
in  each  case  dividing  the  fibres  into  three  equal  lengths  along  the 
circumference.  The  entire  "mat"  could  then  be  removed  from  the  drum  by 
cutting  down  the  centre  of  one  of  the  strips  giving  a  band  approximately 
250  ran  wide  by  450  mm  long  divided  into  three  parts.  This  large  mat  was 
then  cut  into  pieces  57  mm  wide  and  150  mm  long  ready  for  stacking  in  an 
array . 
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3.3.  Moulding  the  composites 

The  mould  used  to  cast  the  specimens  is  shown  in  Fig.l.  This  was  used 
for  unidirectional  specimens  although  other  moulds  were  produced  to 
manufacture  specimens  with  crossplied  reinforcements.  The  end  pieces 
of  the  mould  were  accurately  ground  to  be  57±0‘025  ran  long.  Fibre  mats 
were  cut  to  this  width  on  a  simple  cutter  and  stacked  in  the  mould  with 
brass  or  steel  shims  of  the  correct  thickness  at  each  end  to  space  the 
layers  apart.  After  inserting  the  end  pieces  the  shims  were  strained 
slowly  between  the  grips  of  a  tensile  testing  machine  in  order  to  pull 
the  fibres  tight  and  align  them.  When  full  alignment  was  achieved  the 
end  pieces  were  clamped  tight  and  the  array  was  ready  for  impregnating; 
with  resin. 

The  spacings  achieved  by  this  method  were  not  entirely  perfect.  Figs. 2 
and  3  show  the  best  and  worst  achieved.  However,  compared  with 
commercial  material  (Fig. 4)  there  are  several  advantages.  The  fibres 
retained  their  parallelism  throughout  their  length  and  there  were  very 
few  fibre  to  fibre  contacts,  the  dispersion  being  generally  good.  The 
casts  were  approximately  1  cm  thick.  Twenty  layers  of  fibres  were  used 
in  the  lower  concentrations;  40  or  more  layers  as  necessary  for  higher 
concentrations . 

The  fibre  content  of  Cast  30  was  disappointingly  low  and  so  a  different 
method  was  used  to  prepare  the  remaining  "maximum  concentration"  ca  cs. 
These  were  made  by  packing  a  cylindrical  mould  as  tightly  as  possible 
with  fibres.  The  cast  was  then  infiltrated  with  resin. 

3.4.  Preparation  of  the  epoxy  resin 

The  epoxy  resin  used  to  manufacture  the  specimens  was  Shell  RXE11  with 
Shell  Epikure  RXE11  hardener  which  is  a  mixed  anhydride.  In  order  to 
bring  the  glass  transition  temperature  to  a  lower  value  than  normal, 
for  ease  in  experimentation,  equal  weights  of  an  aliphatic  epoxy  resin 
Araldite  X83/219  were  added  to  the  Shell  resin  before  curing.  In  this 
way  the  glass  transition  temperature  (Tg)  was  lowered  from  about  125°C 
to  65°C  and  the  system  could  be  simply  studied  in  the  laboratory.  Full 
details  of  the  preparation  of  the  various  casts  are  shown  in  Table  1. 

3.5.  Preparation  of  control  specimens 

It  is  possible  that  variation  in  the  chemical  compositions  of  the  resins 
and  changes  in  the  curing  cycle  could  produce  variability  in  the  samples 
which  might  mask  the  effects  of  fibre  reinforcement  alignment  and  fibre 
concentration.  Control  samples  were  therefore  made  from  the  unreinforced 
part  of  each  cast,  being  the  strip  about  15  mm  wide  on  each  side  of  the 
57  mm  central  band.  The  results  obtained  are  shown  graphically  in  Fig. 7 
and  tabulated  in  Table  2. 

3.6.  Design  of  the  test  specimen 

In  order  to  simplify  the  calculation  of  torsional  modulus,  it  was  decided 
to  use  test  pieces  of  circular  cross-section.  The  stress  acting  on  any 
annulus  could  be  assumed  to  be  reasonably  uniform  although  there  would 
of  course  be  variation  in  a  radial  direction.  Test  pieces  were  cut  from 
the  above  casts  at  angles  of  0°,  30°,  45°,  60°  and  90°  to  the  fibre  axis. 
This  enabled  the  moduli  at  various  angles  to  the  fibre  axis  to  be 
determined. 

Test  pieces  were  machined  to  the  shape  shown  in  Fig. 5.  The  splines  at. 
each  end  of  the  specimens  engaged  with  corresponding  pieces  in  the  spring 
loaded  grips  which  are  shown  in  Fig. 6.  These  grips  were  designed  to  hold 
the  specimen  in  shear  while  allowing  for  expansion  and  contraction  during 
the  heating  cycle. 
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3.7.  Measuring  the  torsional  moduli  and  damping  factors 

These  tests  were  performed  with  a  torsion  pendulum  of  conventional  design, 
the  sample  being  immersed  in  a  temperature  controlled  oil  bath.  The 
oscillations  of  the  pendulum  were  recorded  by  means  of  a  suitable  recorder. 

The  frequency  of  oscillation  was  maintained  as  nearly  as  possible  to  1  Hertz 
by  varying  the  inertia  of  the  pendulum  when  necessary. 

From  the  recording  the  torsional  modulus,  'G'  and  the  damping  factor,  "tan  6" 
were  measured.  The  maximum  deflection  of  the  trace  was  never  more  than 
100  mm  at  an  optical  lever  length  of  1  netre.  The  angle  of  torsion  on  a 
28  mm  long  test  piece  was  therefore  only  0*05  radian.  Determinations  were 
made  at  various  sample  temperatures  between  20°C  and  120°C. 

3.8.  Measurement  of  fibre  concentration 

After  the  mechanical  measurements  had  been  made  the  samples  were  cut  up  and 
the  cylindrical  section  used  to  determine  fibre  concentration.  These  were 
calculated  from  the  determinations  of  density  and  by  ashing  the  composite. 

The  results  are  set  out  in  Table  3. 

3.9.  Effect  of  amino-silane  finish  on  the  a  peak 

A  test  to  determine  the  effect  of  the  amino-silane  finish  on  the  position  of 
the  a  peak  was  carried  out.  A  standard  100:100:75  resin/hardener  mixture  was 
cast  in  four  identical  moulds  and  0*8Z,  0*5Z,  0*3Z  and  0*OZ  of  A1100  silane 
added  to  the  four  moulds  respectively.  All  four  were  then  cured  for  16  hours 
at  116°C.  The  test  results  are  plotted  in  Fig. 9. 

3.10.  Effect  of  specimen  length  on  shear  modulus 

To  test  the  effect  of  specimen  length  on  the  measured  values  of  modulus  Cast  17 
was  used.  Specimens  were  cut  with  standard  end  grips  but  with  different  test 
lengths,  namely  55  mm  and  40  mm.  A  40  mm  long  sample  was  also  cut  from  a 
control  specimen.  Shear  modulus  was  measured  at  25°C  and  100°C  and  results  are 
given  in  Table  4. 


4.  DISCUSSION  OF  RESULTS 

4.1.  Fibre  concentration 

An  examination  of  Table  3  indicates  that  the  two  methods  of  determination  of 
fibre  concentration  agreed  to  within  2Z  at  the  higher  concentrations  but  with 
less  accuracy  at  the  very  low  concentrations.  A  comparison  of  Tables  1  and  3 
shows  that  the  fibre  concentration  which  was  achieved  was  generally  fairly 
near  to  the  design  aim. 

4.2.  Results  of  mechanical  tests  on  the  control  samples 
4.2.1.  Variability 

An  examination  of  the  results  in  Table  2  shows  some  variability 
in  the  results  of  the  five  casts  examined.  The  high  values  of 
tan  6  achieved  on  Cast  13  are  probably  attributable  to  the  fact 
that  this  cast  was  submitted  to  a  less  effective  curing  cycle  than 
the  remainder.  The  variability  exhibited  by  the  other  casts  is 
less  easy  to  explain,  however  the  results  indicate  that  Casts  18 
and  20  which  received  the  most  prolonged  curing  cycle  gave  the 
least  variable  results.  Another  source  of  variability  could  be 
the  fact  that  varying  proportions  of  volatile  components  in  the 
resin  boiled  off  during  the  vacuum  treatment.  For  this  reason 
the  procedure  which  was  followed  of  using  as  a  control  for  each 
cast  an  unreinforced  resin  sample  from  the  same  cast,  was  clearly 
correct. 
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A. 2. 2.  Effect  of  the  state  of  cure  on  resin  properties 
Casts  29  to  43  were  poet-cured  at  about  155°C.  This  is  some  20°C 
higher  than  the  earlier  casts.  Comparison  of  the  controls  from 
these  casts  with  those  of  earlier  ones  shows  no  appreciable  change 
in  the  position  of  the  a  peak.  It  is  however  clear  that  the  resin 
mixture  is  probably  not  fully  cured  chemically  and  that  the  addition 
of  more  hardener  would  increase  the  cure  and  thus  the  temperature  at 
which  the  a  transition  occurs.  This  has  an  important  effect  when 
amino-silane  coupling  agents  are  used  as  will  be  discussed  later. 

The  effect  of  increasing  the  amount  of  hardener  is  shown  in  Casts  35 
and  37.  (Pig. 8)  Cast  35  was  prepared  from  proportions  of  100:100:75 
of  the  two  resins  and  hardener  whilst  Cast  37  contained  proportions 
of  100:100:105.  The  o  peak  for  Cast  35  was  at  about  65°C  whilst  that 
for  Cast  37  was  at  76°C.  The  lower  figure  is  within  the  normal 
region  for  this  mixture  when  cured;  the  higher  figure  for  Cast  37  is 
undoubtedly  due  to  increased  cross linking.  It  should  be  noted  that 
both  specimens  were  cured  under  the  same  temperature  conditions. 


A. 3.  Theory 

A. 3.1.  The  elastic  compliance  tensor  for  arbitrary  fibre  direction 
The  symmetry  in  a  unidirectionally  fibre  reinforced  specimen  with 
hexagonal  packing  is  that  of  transverse  isotropy  with  5  elastic 
constants.  Assuming  that  the  1  -  axis  lies  along  the  fibres,  then 
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where  the  S*£j  are  obtained  from  the  S^j  by  rotation  of  axes  1  and  2 

about  the  3  -  axis  (Hearmon  4.4). 
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^  r  2n2(S  -S  )  +  m2S  +  4m2n2(S  -  2S  +  S  ) 


22  23 


55 


11  12  22 


+(m2-n2)2S55  -  [  m3n  (2S12  -  2Sn  +  S^) 


+  mnJ  (2S22  -  2S 


12  ~  SS5>f 


4  2  2  4 

m  S  +  m  n  (2S  +  S  )  +  n  S 

11  12  55  22 


where  m  ■  b 
n  ■  6 


where  0  is  the  angle  between  the  symmetry  axis  and  transformed  axis;  the 
term  in  square  brackets  being  zero  if  free  torsion  is  allowed. 

Therefore,  if  the  elastic  constants  are  known  for  the  unidirectional 
case  then  the  shear  modulus  for  a  sample  cut  at  any  angle  can  be 
predicted. 

In  the  experimental  arrangement  used  the  mounting  of  the  specimen  was 
such  that  bending  was  prevented  and  therefore  the  appropriate  shear 
modulus  is  that  for  the  pure  torsion  (restricted  bending)  case.  GO, 

G30,  G45,  G60  and  G90  have  been  calculated  for  comparison  with 
experimental  results. 

If  further  the  assumption  is  made  that  complex  quantities  can  be 
used  in  the  equations  then  G*  -  Gl  +  GH,  where  Gl  is  the  measured 
modulus  and  G“,  where  Gl  is  the  measured  modulus  and  G'l/G*  »  tan  6 
can  be  calculated  from  the  tensor  of  complex  compliances  S*£j .  The 
validity  of  this  approach  is  discussed  by  Biot  0-1 )  and  Rogers  (12). 

4.3.2.  Calculation  of  the  elastic  compliance  tensor  (Sij) 

The  Sjj  values  have  been  calculated  from  three  different  published 
formulae.  The  simplest  is  the  theory  of  Whitney  and  Riley  (3)  who 
assume  a  composite  unit  cell  of  cylindrical  form  and  apply  classical 
elasticity  theory  to  produce  a  solution  in  closed  form.  This  is 
particularly  convenient  for  calculation  since  it  permits  the  use  of  a 
computer  programme  which  only  requires  input  data  relating  to  the 
elastic  constants  of  fibres  and  matrix  separately  and  the  fibre 
concentration. 

The  second  theory  used  is  that  of  Savin  and  Van  Fo  Fy  (4)  who  give 
approximate  formulae,  derived  from  their  more  detailed  analysis,  for 
the  S{j.  The  third  souce  of  S£j  values  is  from  the  bounds  on  elastic 
moduli  derived  by  Hashin  (13)  ana  by  Hill  (12),  the  validity  of  the 
three  different  approaches  and  their  comparison  with  measured  results 
is  discussed  below. 


4. 3. 2.1.  Whitney-Riley  model 

Whitney  and  Riley  calculate  a  solution  which  is  almost  exact  for  a 
single  composite  cylinder.  (Cylinder  of  matrix  material  concentric 
with  reinforcing  fibre.) 

They  make  one  assumption  in  their  derivation  however,  which  is  that 
the  transverse  Poisson's  ration  v23  is  a  linear  combination  of  the 
Poisson's  ratios  of  fibre  and  matrix  in  proportion  to  fibre 
concentration.  Their  further  assumption  is,  of  course,  that  the 
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single  composite  cylinder  is  representative  of  an  array  of  closely 
packed  cylinders. 

Calculation  of  the  stress  produced  by  stretching  a  uniform  plate 
containing  an  inserted  elastic  disc  shows  that  the  effect  of  the 
disc  will  become  noticeable  at  about  3  radii  *;;om  its  centre.  For 
hexagonal  close  packed  fibres  this  amounts  to  about  10Z  by  volume. 

The  Whitney-Riley  theory,  assuming  that  V23  is  adequately  given  by 
the  linear  approximation  referred  to  above,  should  therefore  be 
valid  up  to  about  10%  by  volume  of  fibre  but  progressively  fail 
thereafter. 

4. 3. 2. 2.  Van  Fo  Fy  and  Savin  formula 

Van  Fo  Fy  and  Savin  set  up  the  equations  of  equilibrium  for  fibres 
in  a  hexagonal  array  using  the  Muskhelishvili  method.  They  chose 
Weierstrass  elliptic  functions,  which  have  the  required  symmetry  of 
the  hexagonal  lattice,  to  construct  a  solution  and  give  the  conditions 
of  compatibility  for  the  matrix  fibre  boundary. 

In  theory  this  method  of  solution  with  the  addition  of  the  solution 
of  the  problem  of  tension-compression  parallel  to  the  fibres  should 
yield  exact  solutions  in  series  for  the  stresses  and  displacements 
in  r.  ,e  body.  They  quote,  however,  only  an  approximate  solution, 
acci  rate  to  15 7.  and  give  the  Sjj  in  terms  of  fibre  and  matrix  moduli 
and  t:he  fibre  concentration.  Their  expressions  for  S55,  and  and 
S  coincide  with  those  given  by  Whitney  and  Riley  but  those  given  for 

and  $2^  differ  considerably.  Since  it  is  these  latter  which 

determine  the  shear  properties  for  off-axis  deformations  the  Van  Fo  Fy 
and  Savin  theory  should  give  better  results  than  Whitney-Riley. 

4. 3. 2. 3.  Hashin  -  Hill  bounds 

Hashin  (13)  gives  bounds  on  the  plane  strain  bulk  modulus,  on  G23  and 
on  Gj2»  while  Hill  gives  bounds  on  v^2  atu*  on  ^1*  ln  hoth  cases  these 
bounds  are  given  for  arbitrary  transverse  phase  geometry  and  therefore 
should  hold  whether  the  packing  is  hexagonal  or  not. 

Values  of  Sjj  for  both  upper  and  lower  bounds  were  calculated  from 
these  formulae  and  used  to  predict  upper  and  lower  bounds  for  the 
shear  moduli  at  angles  to  the  fibres. 

In  all  three  sets  of  formulae  it  is  possible  to  substitute  complex 
quantities  for  the  elastic  moduli  of  fibre  and  matrix  and  therefore 
to  predict  the  modulus  and  damping  of  the  composite.  This  has  been 
done  in  deriving  the  curves.  The  validity  of  this  is  discussed  later. 

4.4.  Comparision  of  theoretical  values  with  test  results 

4.4.1.  Modulus,  glassy  region  (Figs. 10  to  14) 

At  fibre  concentration  up  to  10£,  agreement  between  theory  and 
experiment  appears  to  be  satisfactory.  Both  the  simple  theory  of 
Whitney-Riley  and  the  Van  Fo  Fy  -  Savin  approximation  give  similar 
results  though  differences  start  to  appear  as  the  angle  of  the 
fibres  to  the  specimen  axis  increases  towards  90°  and  brings  in  a 
larger  proportion  of  the  terms  S22  and  The  Hashin  -  Hill 

lower  bound  coincides  everywhere  with  the  Van  Fo  Fy  -  Savin  value. 
Thus,  as  expected  for  fibre  concentrations  where  the  stress  field 
from  one  fibre  cannot  appreciably  affect  the  next  fibre  the 
composite  cylinder  approximation  is  quite  accurate. 
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tor  fibre  concentration*  of  greater  than  10Z  however  the  simple 
theories  do  not  predict  the  shear  modulus  accurately  except  at 
0°  and  90°.  The  Whitney  -  Riley  value  is  in  error  for  the  90° 
case  when  the  concentration  is  over  30Z,  and  in  all  cases  gives 
higher  values  than  the  other  two  theories.  For  the  30°,  45°  and 
60°  specimens  the  measured  modulus  is  greater  than  either  the 
Whitney  -  Riley  or  Van  Fo  Fy  -  Savin  predicted  values  as  soon  as 
the  fibre  concentration  rises  above  10Z  although  the  figures  still 
lie  between  the  Hashin  -  Hill  bounds. 

4.4.2.  Rubbery  region  (Figs. 15  to  19) 

In  this  region  the  differences  between  experiment  and  theory  are 
even  more  marked.  For  all  fibre  concentrations  there  is  poor 
agreement  between  measured  and  calculated  values  of  the  shear 
modulus . 

The  difference  amounts  to  a  factor  of  about  1*5  at  very  low 
concentrations  for  the  0°  orientation  rising  to  3  at  40% 
concentration.  In  all  cases  the  measured  modulus  was  the 
greater  by  this  factor,  which  effectively  rules  out  the 
possibility  of  the  fibres  slipping  in  the  resin,  or  otherwise 
being  ineffective. 

As  the  orientation  of  the  fibres  changes  towards  90°  the 
differences  between  theory  and  experiment  are  reduced  considerably, 
being  noticeable  only  at  the  higher  concentrations  while  there  is 
agreement  at  the  low  concentrations. 

4.4.3.  Damping,  glassy  region  (Figs. 10  to  14) 

For  specimens  with  nearly  rigid  fibres  in  a  more  compliant  matrix 
(Gf/G~60)  it  would  be  expected  that  most  of  the  energy  absorbed 

in  oscillations  would  be  in  the  visco-elastic  matrix  material  and 
hence  that  the  tan  6  curves  for  the  composite  should  always  lie  near 
to  those  for  the  pure  matrix  material.  In  general  this  appears  to 
be  the  case,  the  control  curves  acting  as  good  approximations  to 
the  curves  for  the  reinforced  material  at  all  fibre  angles,  up  to  a 
concentration  of  about  10Z.  It  is  still  true  up  to  40Z  fibre  for 
the  0°  and  90°  cases  but  it  is  very  far  from  being  true  for  the  30°, 
45°  and  60°  cases  at  the  higher  concentrations. 

Theory  would  lead  one  to  expect  a  lowering  of  the  value  of  tan 6  for 
the  off  axis  cases  and  this  is  shown  in  che  curves  drawn  for  the 
various  theories.  The  amount  of  shift  is  however  much  greater  than 
expected. 

4,5,  Possible  reasons  for  the  anomalies 

The  experimental  anomalies  at  30°,  45°  and  60°  for  high  fibre  concentrations 
could  arise  from  several  causes: 

(1)  excitation  of  the  bending  mode  during  testing 

(2)  invalidity  of  the  assumption  that  S12  ■  S21 

(3)  invalidity  of  other  assumptions 

(4)  shift  of  the  a  peak  to  higher  temperatures 

(5)  inapplicability  of  the  Voigt  relation  (Eqn.2) 

These  are  discussed  separately. 


4.5.1.  Excitation  of  the  bending  mode  during  testing 

With  samples  cut  at  angles  to  the  fibre  axis  it  is  difficult  not 
to  excite  the  bending  mode  which  can  naturally  occur  with 
asynmetric  specimens.  In  some  cases  these  oscillations  were 
apparent  but  by  no  means  in  all  and  it  could  not  therefore 
account  for  the  fact  that  at  all  temperatures  in  casts  18  and 
20  t.he  value  of  tan  6  for  the  30°,  45°  and  60°  cases  is  well 
below  the  control  value.  A  full  analysis  of  the  modes  of 
vibration  and  their  interplay  needs  to  be  done,  but  at  first 
sight  one  might  expect  tan  6  to  increase  if  energy  is  being 
absorbed  in  an  additional  mode,  rather  than  to  decrease. 

4.5.2.  Invalidity  of  the  assumption  th.it  Sn  ■  S<n 

It  has  been  assumed  that  S12  -  S21-  Rogers  and  Pipkin  (12>  point 

out  however  that  this  assumption  in  a  time  dependent  system  implicitly 
assumes  the  validity  of  the  Onsager  reciprocal  relations  and  therefore 
needs  experimental  verification. 

For  0°  and  90P  S12  and  S21  occur  as  a  sum  but  for  other  angles  they  do 

not  and  can  play  larger  or  smaller  parts  depending  on  the  angle. 
Hearmon^  notes  a  failure  of  the  relation  in  some  experiments  on 
wood  but  suggests  that  experimental  errors  or  inhomogeneities  may 
be  responsible.  Computed  results  for  S2^  varying  from  0‘5  Sj2  to 

2Sj2  8how  °«ly  changes  in  tan6(Table  5).  It  is  therefore 

concluded  that  this  cannot  be  the  cause  of  the  anomalies. 

4.5.3.  Invalidity  of  other  assumptions 

The  anomalies  in  the  calculated  results  may  also  be  due  to  the 
assumptions  made  in  choosing  input  data  for  the  programmes. 

These  were  (a)  Poisson's  ratio  for  the  resin  ■  0*35 

(b)  The  glass  transition  for  Young's  modulus  (E) 
occurs  at  the  same  temperature  as  it  does  in  shear. 

The  effect  of  the  first  assumption  is  small  as  may  be  seen  from  the 
results  in  Table  6,  calculated  for  Cast  20  at  45°C  and  50°C  using 
E/G  ratios  of  2*7,  2*8  and  2*9  corresponding  to  Poisson's  ratios  of 
0*35,  0*4  and  0*45  respectively. 

G0  is  unaffected  by  changes  in  Poisson's  ratio  since  it  is 
independent  of  it. 

It  is  clear  that  the  changes  in  the  shear  modulus  as  a  result  of 
changing  Poisson's  ratio  are  not  sufficient  to  account  for  the 
anomalies  found. 

It  would  be  of  value  to  have  a  direct  measurement  of  Young's 
modulus,  the  bulk  modulus  or  Poisson's  ratio  at  the  temperature 
and  frequency  of  the  measured  shear  modulus  but  this  is  not  an 
easy  task  if  accuracy  comparable  with  that  possible  in  torsion 
is  desired. 

The  majority  of  methods  (see  e.g.  Neilson  use  forced 

oscillations  or  resonance  at  a  considerably  higher  frequency 
than  that  used  in  this  work  for  shear  measurements.  This 
would  shift  the  relaxation  peak  an  unknown  amount  and  make 
the  use  of  such  results  questionable.  An  attempt  has  been 
made  to  measure  Young's  modulus  and  damping  in  bending  at 
1  Hertz  but  unsatisfactory  results  were  obtained. 
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There  is  some  evidence  (Kono  ^  )  that  the  positions  of  the 

glass/rubber  transition  temperatures  in  shear  and  in  bulk 
differ  but  the  problem  has  been  very  little  studied. 

4.5.4.  Shift  of  q  peak 

Examination  of  the  o-transition  in  its  entirety  shows  that 
there  is  a  lowering  of  the  tan  6  peak  in  the  glassy  region 
accompanied  by  a  general  shift  of  the  transition  to  a 
higher  temperature.  This  could  have  several  causes  but 
the  most  likely  are 

(1)  a  change  of  the  effective  glass  transition 
temperature  by  increased  cross  linking  in  the 
fibre  filled  material 

(2)  a  physical  change  in  this  temperature  by 
reason  of  the  mechanical  interaction  between 
fibres  and  resin 

It  would  seem  that  both  mechanisms  are  playing  a  part  as  will 
now  be  demonstrated. 

4. 5. 4.1.  Cross  linking  of  the  resin 

It  was  mentioned  previously  that  although  the  resin  system 
could  not  be  further  cured  by  more  heat  treatment  it  was  not 
fully- chemical ly  cured.  In  preparing  most  specimens  the  fibre 
was  drawn  through  a  bath  of  10Z  A1100  in  water  and  this  deposited 
a  fairly  large  (but  undetermined)  amount  of  silane  finish  on  the 
fibre.  Sterman  and  Marsden  have  shown  that  the  weight  per 

cent  of  silane  required  is  given  by 

surface  area  of  100  grams  of  filler 

area  covered  by  1  gram  silane 

and  for  the  present  system  this  is  approximately  0*6  weight  per 
cent.  In  section  3.9  it  has  been  shown  that  an  addition  of  this 
amount  of  silane  to  the  same  resin  system  will  give  an  upward 
shift  in  the  glass  transition  temperature  of  4°  to  5°C  and  this 
is  the  order  of  the  shifts  found  in  Casts  18,  20,  30  etc.  which 
use  A1100. 

By  contrast  no  finish  was  used  in  Cast  29  and  the  glass  transition 
temperatures  are  coincident,  within  experimental  error,  with  the 
control  values  (Figs. 20  to  26).  It  seems  clear  therefore  that  the 
shift  in  glass  transition  temperature  between  control  and  filled 
systems  could  be  due  to  the  effect  of  the  silane  finish. 

4. 5. 4. 2.  Interaction  between  fibres  and  resin 

The  results  indicate  another  possible  reason  for  a  shift  of  the 
a  transition.  Casts  33  and  42  were  made  without  silane  yet  a 
clear  shift  of  the  o  transition  relative  to  the  control  samples 
is  evident  in  both  these  casts  (Figs. 27  and  28).  These  two  casts 
contained  higher  fibre  concentrations  and  it  appears  that  this 
shifts  the  o-peak,  broadens  and  flattens  it.  Explanations  for 
this  phenomenon  in  casts  with  a  higher  fibre  concentration  are 
considered  below. 
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4.5.5.  Factors  in  casts  with  higher  fibre  concentrations 
which  could  affect  the  q  transition 

4. 5. 5.1.  Artefacts 

Slippage  of  fibres  in  the  resin  during  the  torsion  test  could 
absorb  energy  and  raise  the  level  of  tan  5.  However  it  is 
difficult  to  see  why  this  should  not  lead  also  to  a  modulus 
much  less  than  the  theoretical  value  instead  of  greater. 

Furthermore  it  should  also  occur  in  the  glassy  region,  yet 
tan  6  is  lower  in  this  region  than  the  control  value. 

4. 5. 5. 2.  Shape  or  type  of  sample 

It  is  not  clear  how  this  could  affect  tan  5  which  is  a  dimensionless 
quantity  although  it  could  conceivably  alter  che  measured  modulus 
values . 

To  test  the  effect  of  specimen  length  on  the  measured  values  of 
modulus  Cast  17  was  used,  specimens  being  cut  with  the  standerd 
end-grips  but  different  test  lengths.  The  results  obtained  are 
tabulated  in  Table  4. 

The  modulus  values  for  the  control  are  quite  comparable  with  those 
for  a  normal  28  mm  specimen  and  the  calculated  moduli  for  the 
glassy  state  (25°C)  agree  reasonably  with  the  measured  values 
although  they  are  a  trifle  low  in  the  case  of  the  long  specimens. 

In  all  cases  however  the  predicted  rubbery  modulus  is  low  as 
compared  with  the  measured  one,  exactly  as  is  the  case  with  the 
normal  specimens. 

4. 5. 5. 3.  Secondary  relaxation  peak 

The  change  in  shape  of  the  a  peak  in  the  highly  reinforced  samples 
could  be  explained  if  there  were  a  high  temperature  absorption 
peak  in  addition  tc  the  normal  one  for  the  flexibilised  resin  at 
65°C.  Such  a  peak  could  be  the  usual  one  for  epoxy  resins  at 
about  120  -  140°C  but  although  it  is  conceivable  that  this  might 
appear  as  a  result  of  cross-linkage  between  resin  and  the  surface 
active  agent  (the  silane  finish  on  the  glass)  it  is  impossible  to 
see  how  this  mechanism  could  exist  when  no  finish  is  present. 

4. 5. 5. 4.  Constraint  of  fibres  on  resin 

Thermal  expansion  could  result  in  compressive  stress  of  up  to  70 
bar  in  the  resin  assuming  perfect  bonding  of  fibres  and  resin. 

The  resultant  restriction  of  free  volume  of  the  resin  could 
increase  Tg  by  several  degrees  (O'Reilly  (18))  but  it  is  unlikely 
that  it  could  account  for  the  very  large  shifts  of  tan  6  found, 
for  example,  in  Casts  34,  42  and  43  (Figs. 27,  28  and  29). 

4. 5. 5. 5.  Fibre-fibre  contact 

We  have  assumed  so  far  the  idealised  model  of  a  uniform  dispersion 
of  fibres.  However  it  is  possible  that  at  high  fibre  concentrations 
appreciable  numbers  of  fibres  may  make  contact  with  each  other  and, 
during  deformation,  cause  a  frictional  contribution  to  the  damping. 

Consider  the  equation  of  motion  of  the  torsional  system 

16  ♦  tee  +  G6  -  0 
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The  solution  for  free  oscillations  is 

6  -  exp  (-tct/2I)  (A  cost  wt  +  B  sin  wt) 

where  -  4  SI  " 

4 12 

The  relations  tan  6  -  X  -  ~  are  used  and  X  measured 

directly  from  the  observed  damping  curve.  There  are  errors 
in  using  this  method  (Struick  (^))  but  they  are  not  large. 

Now  suppose  that  the  frictional  damping  term  <  contains  a 
constant  factor  k 

IC*  +  1c 

then  tan  6  ■  - — - — 

2  Is 

and  the  values  of  tan  6  as  a  function  of  temperature  could 
differ  from  those  expected  theoretically  by  a  constant 
factor.  The  effect  will  be  much  larger  in  the  rubber  region 
where  I  is  small  than  it  is  in  the  glassy  region  where  I  is 
about  300  times  greater.  Comparisons  between  derived  and 
experimental  values  of  tan  6  are  given  in  Table  7. 

The  inertia  terms  in  the  rubbery  region  were  constant  for 
all  three  of  these  casts  so  that  it  is  permissible  to  use 
the  difference  of  tan  6  rather  than  the  difference  of  k. 

Within  the  limits  of  experimental  accuracy  of  tan  6  the 
differences  do  appear  to  be  constant,  which  is  encouraging 
confirmation  of  the  theory  particularly  considering  the 
decade  difference  of  tan  6  between  85°C  and  105°C.  The 
exception  appears  to  be  Cast  43  where  the  differences  become 
greater  at  the  lower  temperatures.  However  if  the  theoretical 
curve  for  tan  6  for  this  case  is  shifted  to  higher  temperatures 
by  the  amount  of  the  shift  found  in  the  glassy  region  then  the 
differences  become  more  regular.  In  considering  the  reason  for 
these  effects  on  tan  5  it  is  supposed  that  it  is  the  friction 
between  fibres  rather  than  fibre  slippage  which  is  important  on 
the  grounds  that  slippage  of  fibres  in  the  resin  should  be  more 
likely  when  no  finish  is  present  than  when  the  fibres  are  bonded 
to  the  resin.  Yet  the  values  of  tan  5  for  Cast  43  at  high 
temperatures  lie  above  those  for  the  (dry)  Cast  42.  Shift  of 
the  curves  for  Cast  42  along  the  temperature  axis  in  fact  give 
reasonable  superposition  on  to  those  for  Cast  43. 

Thus  the  tan  5  curve  in  the  rubbery  region  would  seem  to  be 
related  to  that  for  the  control  resin  (theoretical  curve)  by 
two  factors: 

(a)  The  shift  of  Tg  to  a  higher  temperature  caused 
by  the  chemical  effect  of  the  surface  active 
ager(:  (A1100) 

(b)  The  increase  of  tan  6  by  frictional  work  between 
fibres. 


v  • 
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5.  A  THEORY  OF  FRACTURE  OF  FIBRE  REINFORCED  VISCOELASTIC  MATERIALS 

5.1.  Concepts  of  fracture 

5.1.1.  The  meanings  of  brittle  and  toagh 

Toughness  is  a  remarkably  ill-defined  word  with  many  interpretations 
of  its  meaning.  Intuitively  it  is  felt  that  a  material  which  with¬ 
stands  a  hammer  blow  is  tough;  whereas  one  that  shatters  under  such 
treatment  is  brittle.  This  view,  upon  refinement,  leads  to  the 
various  impact  tests  such  as  Izod  and  Charpy. 

A  second  definition  comts  from  the  load-extension  curve.  If  this 
rises  in  a  more  or  less  straight  line  to  fracture  such  as  in  the 
testing  of  glass,  the  material  is  brittle.  If  there  is  a  yield 
point  and  a  fair  measure  of  extension  after  the  point,  as  in  mild 
steel,  the  material  is  tough.  Under  this  definition  lubber  is  a 
tough  material,  as  are  most  plastics  above  their  glass  transition. 
However,  both  temperature  and  rcte  of  strain  are  very  important  with 
polymeric  materials  and  both  impact  and  stress-strain  methods  can 
give  highly  erroneous  results  unless  this  is  taken  into  account. 

5.1.2.  Modes  of  failure 

In  many  applications  it  is  the  mode  of  failure  which  is  important 
and  for  this  purpose  a  study  of  the  fracture  process  has  proved 
profitable.  This  has  led  to  the  modifications  of  the  Griffith 
crack  theory  for  brittle  solids  to  take  into  account  plastic  and 
other  work  done  in  enlarging  a  crack  in  the  material.  In  metals 
this  extension  was  made  by  Orowan  and  by  Irwin  (22)  and  in 
rubber  by  Rivlin  and  Thomas  (2^),  The  concept  is  of  energy  balance 
between  the  elastic  energy  released  as  a  crack  extends  and  that 
used  in  creating  new  surfaces,  in  plastic  work  or  in  irreversible 
deformation  processes.  If  the  former  exceeds  a  level  known  as  the 
characteristic  surface  energy  or  crack  extension  force  (Irwin)  then 
the  crack  becomes  unstable  and  fracture  occurs.  This  theory  has 
been  highly  successful  in  many  fields  of  materials  and  will  be 
adopted  here. 

5.1.3.  Thermosets 

Thermosets  generally,  are  brittle  rather  than  tough  materials 
although  the  development  of  polyester  and  opoxy  resins  has  given 
industry  materials  which  are  less  brittle  than  the  old  tbermosets, 
but  which  are  still  not  capable  of  absorbing  energy  very  rapidly 
nor  of  sustaining  flaws  under  tensile  loading.  An  indication  of 
the  relative  toughness  of  plastics  and  metals  is  given  in  Tables  7 
and  8. 

It  is  clear  that  in  comparison  with  the  more  common  metals 
thermosetting  resins  do  not  have  very  great  toughness  as  defined 
by  impact  strength  and  surface  energy  determinations  and  even  the 
incorporation  of  reinforcement  does  not  improve  the  figure  by  more 
than  a  factor  of  10. 

5.1.4.  Previous  work 

Various  authors  have  studied  the  effect  of  type  of  hardener,  its 
amount  and  the  temperature  of  curing  on  the  mechanical  properties 
of  epoxy  reains.  Lewis  and  Ramsey  using  Epon  828  and  diethylene 
triamine  found  fracture  energy  (area  under  stress-strain  curve) 
increased  with  curing  by  a  factor  of  about  4.  This  is  in  accordance 
with  common  experience  that  a  room  temperature  amine-cured  epoxy  is 
much  more  brittle  than  one  cured  at  a  higher  temperature.  Mostovoy 
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and  Ripling  (27)  confirm  this  using  a  measure  of 
*  Gc  ■  (Pc2/2b)|£  where  "PJJ  is  the  applied  load  at  the  instant 

of  fast  crack  extension,  "b"  specimen  thickness,  "C"  specimen 
compliance,  and  "a"  the  crack  length.  They  used  DER  332  and 
'  TEPA  and  showed  that  G0,  the  characteristic  surface  energy 

increased  with  increased  post  cure  temperature  and  also  with 
increasing  amounts  of  TEPA  until  a  limit  was  reached. 

Burhans  et  al(28)  have  shown  that  decreasing  the  distance 
between  "active  sites"  (the  carbon  atoms  furthest  from  the 
ether  linkage  in  epoxies)  increases  both  modulus  and  strength 
but  this  is  generally  at  the  expense  of  toughness.  Conversely, 
toughened  epoxies  can  be  made  by  incorporating  long  chain 
molecules  of  either  resin  or  hardener  into  the  cross  linked 
network.  Thus  the  long  chain  aliphatic  epoxy  Araldite  X83/219 
imparts  toughness  to  the  systems  in  which  it  is  mixed  and  the 
same  effect  is  found  by  using  dodecynyl  succinic  acid  as  hardener 
with  a  standard  bisphenol  A  or  epoxy  novolac  resin  system.  There 
is  considerable  evidence  that  the  curing  of  epoxy  resins  may  be 
excessive  if  maximum  toughness  is  required.  Lewis  and  Ramsey 
found  that  cure  at  60°C  for  four  hours  followed  by  a  short  post¬ 
cure  at  190°C  gave  the  best  results,  and  it  is  probable  that 
this  left  an  appreciable  number  of  groups  unreacted.  It  can  be 
seen  that  resins  are  no  exception  to  the  general  rule  found  in 
metals  that  high  strength  gives  low  toughness  (or  in  other  words 
poor  resistance  to  crack  propagation)  and  vice  versa.  Thus  the 
strongest  steels  are  the  most  brittle  and  the  same  is  true  of 
resins.  (Vincent  (36) t  Bobalek  and  Evans  (37)), 

5.1.5.  Effect  of  fibre  reinforcement 

Fibre  reinforcement  of  thermosetting  resins  is  claimed  to 
combine  strength  and  toughness.  In  order  to  investigate 
this  claim  the  mode  of  failure  of  a  fibre  reinforced  resin 
must  be  considered  in  some  detail  taking  into  account  the  way 
in  which  elastic  strain  energy  can  propagate  through  the  material 
and  in  what  way  it  can  be  dissipated. 

5,2,  Theoretical  considerations  of  fracture  ir.  a  fibre-reinforced 
resin  under  tension 

As  a  tensile  specimen  with  unidirectionally  oriented  fibres  is  loaded 
elastic  energy  is  stored  in  the  fibres  and  the  resin  while  a  certain 
amount  is  also  dissipated  in  the  resin  (since  it  is  viscoelastic).  At 
a  certain  strain  either  the  resin  will  fail  from  some  stress  concentration 
such  as  a  bubble  or  one  fibre  will  break.  The  result  will  be  similar  in 
the  two  cases  in  that 

(a)  the  remainder  of  the  fibres  will  carry  a  slightly  increased 
load 

and 

(b)  a  stress  concentration  will  exist  at  the  tip  of  the  crack 
and  will  therefore  increase  the  likelihood  of  its  spreading. 

5.2.1.  Energy  at  fracture  site 

Three  conditions  are  possible!  either  the  crack  will  spread 
under  the  applied  load,  or  a  new  crack  will  develop,  or  the 
system  will  still  be  stable.  To  ascertain  which  will  occur 
it  is  necessary  to  calculate  whether  the  strain  energy  released 
by  spreading  of  the  crack  (including  the  breaking  of  fibred  is 
sufficient  to  supply  the  needs  of  the  new  surfaces,  and  of 
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viscoelastic  and  frictional  work.  The  fact  that  the  fibres  may 
not  have  a  uniform  strength  or  breaking  strain,  but  a  variable 
one  distributed  about  a  mean  value  must  also  be  taken  into  account. 

In  addition  it  is  unlikely  that  all  fibres  are  equally  well 
aligned  or  that  they  are  equally  strained;  th’*s  will  also  produce 
variability. 

5.2.2.  The  proposed  model 

In  setting  up  the  proposed  model  three  assumptions  have  been  made. 

5. 2. 2.1.  Assumption  1 

It  is  assumed  that  when  one  fibre  breaks  the  disturbance  in  local 
stress  extends  to  nearest  neighbours  only  and  for  a  length  on  either 
side  of  the  break  equal  to  the  elastic  transfer  length  (defined  below). 

5.2.2. 2.  Assumption  2 

It  is  also  assumed  following  Kelly  (29)  and  Cox  (30),  that  the  shear 
stress  in  the  resin  at  the  surface  of  the  broken  fibre,  at  a  distance 
x  from  the  broken  end  is 


t(x) 


2*n  R/rf 


e  sinh  3  (V2  ~  x)  .  (i) 

cosh  3  *72 


where  is  Young's  modulus  for  the  fibre 

Ggj  is  shear  modulus  for  the  resin 
e  overall  strain  in  the  composite 
2R  mean  separation  of  fibres  laterally 
rj  fibre  radius 
1/ 2  elastic  transfer  length 
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Together  with  this  expression  for  the  shear  stress  in  the  matrix  the 
stress  in  the  broken  fibre  at  a  distance  x  from  the  free  end  is  taken 
to  be 


abf 


Ef  e 


[* 


cosh  3  (V2  “  x) 
cosh  3  V2 


■] 


(ii) 


for  o<x<  *7 2 


It  it  clear  that  this  expression  can  uever  equal  the  stress  in 
the  unbroken  fibres  unless  i  is  very  large. 

However,  if  0*/2  ■  2  the  stress  at  x  ■  *72  will  reach  70*  of  the 
unbroken  stress,  if  0*72  -  3  it  reaches  90*.  We  therefore  define 
*72  the  elastic  transfer  length  to  be  equal  to  3/e  .  (iii) 

5. 2. 2. 3.  Assumption  3 

The  shear  stress  in  the  resin  at  point  (r,x)  is  taken  to  be: 


t  (r,x)  - 


rf  t  (rf  ,  x) 


/  Ef  Gm  e  rf  sinh  0  (£/2  -  x) 
2  £n  R/r r  r  C08h  6  lh 


Hence  the  load  in  an  adjacent  fibre,  P,  is  related  to  this  shear 
stress  by: 

£  ■  2,rt Tf 

Consider  the  stress  aaf  ■  -JL.  at  an  adjacent  fibre  distant  2R 


from  the  broken  one.  Then: 


7?  2’  rf  | 


T  (2R,x)dx 


'Ef  &m 


e  rf  cosh  0  (*/2-x) 


f  /  2£n  R/rf  C08h  0  lh 


The  constant  C  is  found  by  equating  <>af  to  0  at  x  •  *7 2  since  we 
are  only  interested  in  the  incremental  stress 


0  - 


C+.L  /Ei-9a _ i - 

PR  / 2*n  R/rf  cosh  0  */ 2 


0 R cosh  0  */2  /2*n  R/r 


G”—  [  cosh  0 


(‘/rx)-l] 


rf  /  Ef  £n  R/r. 
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But,  by  definition,  W/2  ■  3,  hence  cosh  &V2  “  10 

**f(o)  -  (°*45  rf/R)  .  (V) 

The  energy  released  when  a  fibre  breaks  can  now  be  calculated  and  it 
can  be  determined  whether  the  stress  on  an  adjacent  fibre  will  be 
sufficient  to  break  it. 

5.2.3.  Energy  changes  during  fracture 

Only  the  "transfer  zone"  mentioned  above  is  being  considered,  i.e. 
the  energy  available  for  spreading  of  the  crack  comes  only  from  a 
re-arrangement  of  the  stresses  in  the  transfer  zone  and  not  from 
any  relaxation  of  stress  cutside.  The  elastic  strain  energy  in  the 
zone  before  fibre  fracture  is  '$  <rc  cir (R) ■  ec.  After  fracture  it 
is  the  sum  of  Ebf  (the  strain  energy  in  the  two  broken  portions  of 
the  fibres)  +  Em  (the  energy  in  the  matrix  material) ,  and  the  energy 
released  on  fibre  fracture  is  therefore  Ec  -  Ebf  -  Em-  Expressions 
for  Ebf  and  Eq  are  developed  below. 

It  is  assumed  that  this  energy  goes 

(1)  to  increase  the  strain  energy  of  the  six  adjacent  fibres 

(2)  to  supply  the  surface  energy  of  the  new  surfaces  created 

(3)  to  supply  the  heat  and  other  energy  absorbed  during  fracture 


5.2.4.  Energy  expressions 

5. 2. 4.1.  Energy  in  the  broken  fibre 
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lit2  E  2  e2  j  _  co»h  B  (&/2-x) 
Ef  *  cosh  6  */2 


r2e2Ef  4/2  -  2/g  tsnh 


s«/2*(‘/4.l^J-i).«ch26f]...  (»i) 


5. 2. 4. 4.  Energy  in  the  matrix  material 


Em-2  *J 


/R  /V2  2 
/  /  LiH>  2  ir  r  dr  dx 

rf  4 


IIr2/  / sinh2  6  (*/2-x)  dx 


'f  ° 


where  K 


.  LZ. LJS - 1 .  — 

h  *n  R/rf  cosh  B  */2 


Em-  *n  R/r  J  J 2  (cosh  2g(4/2-x)  -i) 

4 

■  £  k2  l”  S  [-isii-e-A  -  v4  ] 

ie  .  .  *  Ef  rf2  e2  fsinlLM.  -  */.  1 

^  cosh2  ji.  L  4  B  4  J 


(vii) 


5.2.4. 5.  Energy  in  adjacent  fibres 

!lh  . 


E#£  ■  6  x  2 it  rf2  *  !  /  dx 


6  tt  rf2  of2  /  t^2  [cosh  g  (*/2-x)-lI 

Ej  /  4R2  cosh 
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mm 


.  6  1  2  qf2  £12  sech2  0  A/2  f31^  +  “SiLJLl  -  2.  “S.h ■£ f/J  (viii) 

Ef  4R2  L  4  g  g  J 


5. 2. 4. 6.  Energy  released 
Energy  released  -  Ec  -  E^-  Effl 


-  oc  tv  R2  £  -  ir  rf2  e2  Ef  £*/2  -  2/g  tanh  0  */2  ♦  -”h  g  -  «ech2 

-  i  e2  Ef  n  rf2  £  ♦  \  e2  £m  ,  <R2-rf2)  £ 

2  ^  Ef  [*72  “  7  tanh  ^2  ] 


-  «  r< 


But  0  -  |-  tanh  ^  =| 


Energy 


released  -  n  e2  £  rf2  Ef  -  ♦  J  E,  ^  -  J  Emj 


2  ,  2 
ir  e  £  rf 


!u 


* 


R2-rf2 

'f2J 


(ix) 


I 


5.2.4. 7.  Rate  of  release  of  energy 

Normally  Ef  >  10  En  so  the  first  term  predominates. 

*  2 
Hence  in  the  Griffith  -  Irwin  sense,  for  extension  of  the  crack  by  area  ir  R 

the  rate  of  release  of  energy 


Gc  -  55. 

c  6A 


(x) 


Note  that  this  is  more  than  l/3  of  the  originally  stored  energy. 


5. 2. 4. 8.  Dissipation  of  released  energy 
The  energy  released  in  equation  ix  is  used 

(a)  to  increase  the  strain  in  adjacent  fibres 

(b)  to  suoply  surface  energy  for  the  new  surfaces  created 

(c)  to  supply  heat  and  acoustic  energy  and  possibly  energy  for  molecular 
deformation  processes  such  as  orientation. 
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The  amounts  are: 


Eaf  " 


ii  it‘  E  ,eCh2  £J  fi*  *  ».i«b  e »  .  2  »i°h  8  vil 

2  «2  1  2  [4  4  6  8  J 


2l  £f4  -2  F  » 
80  R?  e  Ef  4 


#  o 

Surface  energy  for  fibres  •  2  y^i  rf 
Surface  energy  for  resin  ■  2  Yffl1T  (R2  -  rf2) 


Dissipated  energy  ■ 

_  y  p 

(Yf  *  560  ergs  /cm2  for  glass.  ym  *  2  x  10^  to  2  x  103  ergs/ctn  for 

epoxy  resin).  Thus  the  energy  equation  becomes  Energy  released  -  Energy 
absorbed,  i.e. 


"  [lo  ^  *2  Ef  1  *  2^f  **f2  ♦  2Tm  "  (R2  -  rf2)  ♦  Ed 
Dividing  by  it  R^  £  gives  the  energy  per  unit  volume  of  the  transfer  zone: 


[?■ • 1  ■■  } [i (# •* - 

•W'-fifi)  -ftr] 


5.2.5.  Criterion  for  fracture  of  neighbouring  fibre 


It  is  now  possible  to  decide  whether  this  difference  is  sufficient 
to  break  a  neighbouring  fibre. 

For  the  six  nearest  neighbours  of  the  broken  fibre  the  probability 
that  under  a  strain  e  five  will  remain  unbroken  and  the  sixth  will 
break  when  a  strain  increment  de  is  applied  is: 


6  [l  -  F  (y)]  f  (y)  dy 


where  y  *  F(y)  the  cumulative  distribution  function, 
o 

f(y)  ■  —  and  a  the  standard  deviation, 
dy 
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P 


Hence  the  expectation  of  another  broken  fibre  in  this  set  of 
six  will  be 

36  [l  -  F(y)]  f(y)  dy 


and  equating  this  to  unity  we  find  the  condition  that 

dy  - - - - r - 

36  [1  -  F(y)]5  f(y)  . 

for  one  of  the  six  to  break. 


(xi) 


Example  Now  dy  ■  ^ 


6  ['  - F  (Hr)]5  '  (Hr)  f  • 1  <»“> 

Table  9  gives  calculated  values  of  61  l-F(y)]  5  f(y)  for  values 
of  y  from  -  »  to  +  1. 

An  examination  of  the  values  shows  that  the  maximum  is  at 
about  y  *  “1  and  is  0*611. 

1 


Hence  dy  = 


-  0*273 


6x0*611 

at  y  -  -1  e-oy  +  ?  •  e  -  o 


and  de  =  £  ■  a  x  0*273 
6 

or  r.  ■  1*638  0 
i.e.  c  ■  2*638  a 


Hence  if  the  coefficient  of  variation  is  less  than  38Z  and 
the  fibre  is  strained  to  within  one  standard  deviation  of 
the  mean,  catastrophic  failure  may  occur. 

If  y  ■  -2,  c  ■  5*47  a  for  catastrophic  failure 

That  is,  if  the  coefficient  of  variation  is  lower  than  18*3Z 
the  strain  must  not  exceed  two  standard  deviations  from  the 
mean. 

5.2.6.  Condition  for  crack  propagation 

The  condition  dy  - - - - *  .  ~  for  failure  of  one  of 

36  [l-F(y)P  f(y) 

the  six  adjacent  fibres  implies  a  value  of  the  incremental 
energy  which  can  then  be  used  as  a  criterion. 

In  a  system  with  strain  energy  J  Ec^  per  unit  volume  an 
increase  of  strain  to  e  +  de  means  an  increase  in  strain 
energy  density  by  Eede.  Hence,  if  y  is  such  as  to  satisfy 
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Che  equation  xi  above  then,  since  y  *  Ini  we  have  Che  condition 

a 

for  propagation  of  a  crack  that  the  energy  released  on  breakage 
of  one  fibre  shall  exceed  the  incremental  strain  energy  derived 
from  the  term  Eftde  at  some  point  on  a  fibre.  At  this  stage  we 

shall  assume  that  it  is  not  the  total  energy  in  the  transfer  zone 
but  the  energy  density  at  a  point  on  an  adjacent  fibre  that  is 
the  critical  factor.  Hence  our  criterion  for  failure  is 

Energy  released  *  Energy  absorbed  >  Eje6e  .  (xiii) 

where  e  and  de  are  related  by  the  statistical  condition  above, 
that  is, 


6e 


36  [l-F  (*J)]»  f  (M) 


5.2.7.  Time  dependence  of  modes  of  failure 
5. 2. 7.1.  Elementary  example 

In  a  bundle  of  NQ  fibres  the  expectation  of  failure  up  to  and 
including  y  is  N  -  N0  F(y)  where  y  -  £p- 


dN 


*»  5,#ofWj|f 


•  N„  f(y)  i  if  a£ 


a 
vdt 


dN  -  n0  f(y)  0 


Hence  the  average  time  for  one  fibre  to  break  will  be 


dt 


N_  f(y)  v 


(xiv) 


If  this  time  exceeds  the  retardation  time  T  of  the  resin  the 
stress  in  the  resin  will  have  time  to  relax,  reducing  the 
strain  energy  available.  If  it  is  shorter  than  T  then  strain 
cannot  relax  before  a  fibre  is  due  to  break,  and  strain 
magnification  will  be  present. 

For  the  six  nearest  neighbours 
AN  -  36  [l  -F(y)]5  f(y) 

and  |j:  -  36  [l  -F(y)]  f(y) 


Ay 


v 

a 


Hence  for  the  breakage  of  one  fibre  of  the  six,  a  time  At 
ia  required 


36  [1  -  F(y)]5  f (y) 


In  this  time  the  energy  density  will  change  as  the  modulus 
of  the  resin  changes  from  its  unrelaxed  to  a  more  relaxed 
state  according  to  the  relation  (for  a  simple  spring  dashpot 
model) : 

_At_ 

E(At)  -  ER  +  (Eu  -  ER)e  “'T 


Hence  the  inequality  (xiii)  becomes 


This  inequality  gives  the  strain  for  brittle  failure  in  terms 
of  parameters  of  the  composite. 

Clearly  brittle  failure  is  more  likely  when  Ef  is  large,  when 
T  is  small,  or  when  v  or  T  are  large.  The  equation  also  takes 
into  account  the  proximity  of  e  to  e  the  mean  breaking  strain 
of  the  fibres  through  the  presence  of  f(y)  in  the  denominator 
of  the  right  hand  side.  The  retardation  time  T  is  of  course 
temperature  dependent,  low  temperatures  giving  long 
retardation  times. 

For  maximum  toughness  a  medium  to  low  proportion  of  fibres, 
preferably  of  very  variable  strength,  and  of  large  diameter 

(since  this  increase  V2  the  transfer  length)  is  required. 
Supporting  evidence  for  this  comes  again  from  work  on  fibre 
reinforced  metals  by  Cooper  and  Kelly  (42). 

The  importance  of  the  rate  of  testing  is  shown  by  McAbee  and 
Chmura  (40)  who  tested  fibreglass  laminates  at  varying  speeds, 
finding  higher  strength  and  modulus  at  the  higher  rates  of 
strain. 

The  author  (50)  has  shown  that  the  mode  of  failure  of  a 
longitudinally  reinforced  resin  sample  when  tested  in  tension 
depends  critically  upon  the  rate  of  testing;  slower  rates  led 
to  a  more  fibrous  type  of  fracture  while  high  speeds  gave  a 
brittle  type. 
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An«  effect  of  fibree  upon  each  other  during  fracture  has  been 
eeen  to  be  the  same  whether  in  a  fibre  reinforced  material  or, 
by  considering  molecular  chains  as  "fibres",  in  a  r*,*jJ* 
same  effect  is  found  in  the  tearing  of  textiles  in  work  publi«hed 
by  Taylor  ,  O'Brien  and  Weiner  (44),  Spencer-Smith  (45)  and 

Kilbv  (4b).  The  general  conclusion  is  that  whon  interfibre 
effects  are  small  (low  twist,  low  interfibre  fiction)  the  tear 
strength  is  high  while  if  the  fibres  are  highly  twisted  or  woven 
in  such  a  way  that  friction  is  high,  or  coated  with  crease 
resistant  or  fire  resistant  finishes,  then  the  tear  resistance 
will  be  low.  This  latter  case  corresponds  to  the  fibre  reinforce¬ 
ment  condition  of  a  short  transfer  length. 

It  is  not  yet  possible  to  answer  qualitatively  the  question  "Are 
strength  and  toughness  incompatible?"  To  do  this  a  more  careful 
study  of  the  energy  changes  in  fracture  of  a  composite  is  needed, 
taking  into  account  all  energy  absorbing  mechanisms  such  as  the 
one  proposed  by  Cook  and  Gordan  (  '  ■ 


5,3.  Predictions  from  the  theory 

Using  equation  (xvi)  as  a  criterion  of  brittle  fracture,  calculations 


resin  system.  The  input 

data 

were: 

Ef  -  7*3  x  1010 

N/m2 

Gr  (resin)  ■  2*1  x  106 

N/m2 

Gy  (resin)  «  l«05x  109 

N/m^ 

yf  (glass)  »  0-56 

J/m^ 

Ym  (resin)  •  2  *  10 

J/m^ 

Concentration  of  fibres: 

10%  to  70%  by  volume 

Strain  rates: 

10"5,  10'3,  10"1 

Retardation  times: 

1000,  1,  1/ 1000s 

Mean  fibre  breaking  strains  e: 

0.0001,  0-01,  0*1 

Standard  deviations: 

0*le,  0*5e 

The  strain  e  was  then  varied  from  zero  to  e  in  steps  of  0*1  e  and  by 
linear  interpolation  the  value  of  e  for  catastrophic  failure  was  then 

found. 

The  results  are  given  graphically  in  figures  31  to  34.  Examination  of 
these  figures  indicates  the  following  conclusions: 

(a)  Catastrophic  failure  should  never  occur  when  T  -  0*5  e 

(b)  Catastrophic  failure  should  never  occur  when  c  -  0*001  (0*1%  strain) 

(c)  For  e  -  0*01  (13!  strain)  (Fig. 31),  the  strain  at  which  catastrophic 

failure  will  occur  falls  as  the  concentration  of  fibres  is 

increased  from  about  0*9  l  at  20%  concentration  to  0*81  e  at  70%. 
Below  20%  concentration  catastrophic  failure  should  not  occur. 
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The  curve  for  short  time-constant  (resin  in  rubbery  state)  lies 
below  the  curve  for  longer  time-constants.  Thus  if  the  resin  is 
in  the  glassy  state  or  a  rapid  strain  rate  is  employed  this 
results  in  a  higher  breaking  strain  than  occurs  in  the  rubbery 
state. 

If  the  size  of  the  reinforcing  fibre  is  10  micron  (Fig. 32) 
instead  of  100  micron  theory  indicates  that  catastrophic  failure 
should  not  occur  except  for  the  rubbery  state  (T  ■  1/1000)  and 
for  the  slowest  strain  rate  in  the  intermediate  state  (T  ■  1). 
Thus  on  this  theory,  for  medium  strength  fibres  (c  ■  0*01)  the 
smaller  fibres  result  in  a  tougher  material. 

(d)  Fore  -  0*1  (10Z  strain)  (Fig. 33),  the  fibres  are  of  maximum 
strength  and  the  pattern  of  behaviour  is  somewhat  different 
from  that  found  for  low  strength  fibres.  The  same  trend  to 
lower  breaking  strains  at  higher  fibre  concentrations  is, 
however,  found  but  the  curves  for  the  rubbery  state  now  lie 
above  those  for  the  glassy  state  and  the  curve  for  the  lowest 
strain  rate  lies  above  those  for  the  higher  rates,  for  any  one 
choice  of  T.  Thus  a  material  reinforced  with  very  strong, 
uniform  fibres  (o  -  0*12)  will  be  tougher  at  low  rates  of  strain 
and  in  the  rubbery  state  (vT  low)  than  at  high  rates  of  strain 
or  in  the  glassy  state  (v  T  high).  There  is  little  change  with 
the  radius  of  the  reinforcing  fibres,  the  curves  for  rf  -  10y 
(Fig. 34)  being  very  similar  to  those  for  r^  -  100. 


5.4.  Discussion 

The  foregoing  theory  is  an  attempt  to  introduce  quantitative  relationships 
into  the  subject  of  the  fracture  of  fibre  reinforced  resins  and  to  take 
into  account  rate  effects. 

It  has  already  been  noted  by  Kelly  (38>  that  too  high  a  concentration  of 
fibres  in  a  metal  matrix  gives  a  brittle  composite  and  the  same  is  true 
in  a  resin  reinforced  by  fibres,  (Wadsworth  and  Spilling  (51).  The  analogy 
with  the  effect  of  cross  linking  on  the  strength  of  rubbers,  analysed  for 
example  by  Bueche  (39),  is  close  and  it  is  suggested  that  this  is  because 
the  theory  required  is  the  same  in  the  two  cases. 

Experimental  work  shows  that  a  highly  cross  linked  resin  is  more  brittle 
than  one  in  which  a  proportion  of  possible  cross  links  have  not  been  formed. 
Thus  Lewis  and  Ramsey's  (26)  results  and  the  findings  of  Van  den  Boogaart  (41) 
together  with  the  work  of  Burhans  (28)  suggest  that  maximum  toughening  in  a 
pure  resin  occurs  when  it  is  undercured.  Can  these  results  be  explained  by 
considering  a  polymer  network  as  though  it  were  a  molecular  scale  fibre- 
reinforced  material?  This  is  effectively  the  basis  of  Bueche 's  treatment 
of  a  cross  linked  network  and  it  is  obviously  qualitatively  satisfactory. 

It  is  proposed  to  add  to  Bueche' s  theory 

(a)  the  concept  of  transfer  length 

(b)  the  consideration  that  chains  may  have  very  variable  breaking 
extensions 

(c)  the  retardation  time  of  the  resin. 

From  the  latter  concepts  it  is  predicted,  using  calculated  results  from  the 
reinforcement  theory,  that  if  the  breaking  extensions  of  the  chains  are 
very  uniform  brittle  failure  is  to  be  expected  (glassy  stste)  whereas  if 
they  are  variable  it  is  expected  that  the  failure  of  one  chain  will  be 
accommodated  by  its  neighbours  in  a  reasonable  time. 
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The  time  to  be  considered  is  that  for  the  relaxation  of  stress  in  the 
neighbourhood  of  a  broken  chain.  The  relevant  relaxation  time  will 
therefore  be  that  for  a  small  chain  segment,  that  iR  the  y  relaxation 
at  low  temperatures.  Neilsot  (31)  states  that  at  temperatures  well 
below  Tg  two  rigid  polymers  can  differ  greatly  in  impact  strength. 

For  example,  cellulose  nitrate  and  polycarbonates  show  much  higher 
impact  strengths  than  polystyrene  or  PMMA:  Neilsen  holds  that  this  is 
due  to  the  secondary  transitions  in  the  former  polymers  (see  also 
Bobalek  and  Evans  (48)). 

In  epoxy  resins  there  are  two  well  defined  relaxation  regions.  The 
crrelaxation  or  glass  rubber  transition,  depending  upon  the  hardener 
used  and  the  resin  system,  lies  above  100°C  and  is  the  main  transition 
region.  The  Y  or  low  temperature  transition  lies  below  -50°C  and  is 
associated  with  the  opening  of  the  epoxy  ring  to  form  a  bond  with  the 
hardener. 

The  curing  of  an  epoxy  increases  the  intensity  of  this  relaxation  and 
shifts  it  (in  an  amine  cured  resin)  from  about  -75  C  to  -50  C.  The 
use  of  a  flexible  hardener  (e.g.  dodecynyl  succinic  anhydride)  moves 
the  y  peak  to  -90°C  and  introduces  a  subsidiary  6  peak  at  -160°C 
(Speake  (49)).  The  same  shift  of  the  Y  peak  to  lower  temperatures  is 
found  by  the  incorporation  of  Araldite  X83/219  (diglycidyl  ether  of 
polypropylene  glycol)  into  the  network  structure  and  in  this  case  also  a 
6  peak  starts  to  appear  at  about  -140°C.  Chemical  undercure  has  also 
been  shown  to  move  the  Y  peak  to  lower  temperatures.  Undercured 
materials  are  tougher  than  resins  made  without  the  flexible  groups, 
the  impact  strength  being  increased  by  a  factor  of  about  2. 

We  may  be  able  to  explain  the  increased  toughness  (or  resistance  to 
catastrophic  failure)  of  these  resins  by  using  the  fibre  reinforcement 
model  and  taking  the  relaxation  time  relevant  to  the  y  relaxation.  This 
is  a  low  energy  transition  and  moves  about  6  to  10°C  per  decade  of 
frequency.  The  results  quoted  were  obtained  at  1  Hertz  so  that  at  room 
temperature  the  -50°C  relaxation  should  appear  at  about  1  megahertz. 

Thus  the  relaxation  time  at  room  temperature  can  be  taken  as  10-6  second. 
For  the  lower  temperature  relaxations  found  when  flexible  groups  are  used 
in  either  resin  or  hardener  the  relaxation  time  may  well  be  io-9  second. 

Thus  at  reasonable  strain  rates  the  product  vT  in  equation  (xvi)  will  be 
very  small  and  the  material  will  behave  in  a  rubbery  fashion.  That  is, 
the  breakage  of  one  chain  will  be  rapidly  accommodated  by  the  neighbouring 
chains  and  the  transfer  zone  will  be  large,  the  crack  blunted.  At  higher 
rates  of  strain  vT  may  approach  glassy  values  if  T  is  only  io“6  but  is 
unlikely  to  if  it  is  10“9.  Thus  under  impact  conditions  the  material  is 
likely  to  be  tougher  when  the  y  peak  lies  below  50°C. 

It  may  be  possible  to  be  more  quantitative  in  the  use  of  the  fibre 
reinforcement  analogue  and  to  design  a  resin  for  toughness  and  strength 
possibly  by  choosing  copolymer  blends  to  give  the  right  balance  between 
load  bearing  chains  and  energy  absorbing  ones. 

One  may  expect  the  pure  resin,  pure  hardener  system  to  be  virtually 
useless  because  of  brittleness  but  that  a  distribution  of  molecular 
weights  and  a  mixture  of  different  hardeners  would  give  optimum  properties. 
Much  remains  to  be  done,  however,  before  such  molecular  engineering  is 
possible. 


CONCLUSIONS 


6.1.  Modulus  in  the  glassy  region 

Within  the  glassy  region  the  agreement  between  predicted  and  measured 
values  of  the  shear  modulus  is  good.  Since  predicted  values  for  specimens 
cut  at  various  angles  to  the  fibres  axis  involve  all  five  predicted 
elastic  constants  it  can  be  inferred  that  all  five  are  reasonably  accurate. 
The  experimental  values  all  tend  to  lie  above  the  predicted  ones,  however, 
in  agreement  with  the  comment  by  Chassis  and  Sendeckyj (8) . 

6.2.  Modulus  in  the  rubbery  region 

In  the  rubbery  region,  good  agreement  is  not  obtained.  It  is  difficult  to 
account  for  this  discrepancy  other  than  by  assuming  that  it  is  the  basic 
assumption  of  the  theories  which  are  in  error.  Thus  it  is  by  no  means 
certain  that  reain  in  the  interstices  between  fibres  has  the  same  elastic 
properties  that  it  has  in  che  bulk.  There  are  also  unknown  compressive 
stresses  involved  which  may  play  a  part.  However,  it  is  likely  that  the 
largest  error  consists  in  using  a  homogeneous  anisotropic  model  for  a 
heterogeneous  material  and  it  is  possible  that  a  much  larger  specimen 
would  have  approximated  better  to  the  homogeneous  ideal. 

6.3.  Tan  6 

Turning  to  the  loss  factor  the  epperiments  show  clearly  that  for  the  system 
used  it  is  perfectly  valid  to  substitute  complex  quantities  for  the  elastic 
moduli  in  explicit  formulae  such  as  those  of  Whitney-Riley,  Hashin  or 
Van  Fo  Fy  and  Savin,  and  that  the  predicted  tan  5  curve  corresponds  closely 
to  that  observed  provided  that  chemical  interaction  of  surface  active  agents 
with  the  resin  and  secondary  sources  of  energy  absorption  are  taken  into 
account. 

A  comparison  of  our  results  with  predictions  based  on  other  theories  now 
available  may  make  it  possible  to  state  which  explicit  formulae  are  the 
best  to  use  for  predicting  the  complex  moduli  of  a  fibre  reinforced  system. 
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TABLE  1 

DETAILS  OF  CASTS 


CAST  No. 


13 


*  15 


Unknown 


NOMINAL 

CONTRACTION,  v/v 

FINISH 

APPLIED 

32 

Yes 

32 

m 

122 

Yea 

122 

Yes 

122 

Yes 

402 

Yes 

402 

No 

Maximum 

Yes 

Maximum 

No 

Maximum 

Yes 

0 

No 

0 

No 

Maximum 

No 

Maximum 


Control  for  42  and  43 


CURE 


24  hrs  130°C 


16  hrs  100°C 
►24  hrs  140°C 


16  hrs  100°C 
►24  hrs  130°C 


17  hrs  105°C 
►24  hrs  130°C 


24  hrs  110°C 
►24  hrs  130°C 


20  hrs  110°C 
24  hrs  130°C 


17 i  hrs  98°C 
►  9  hrs  115°C 
►24  hrs  145°C 


16  hrs  110°C 
►26  hrs  155°C 


17  hrs  120°C 
►26  hrs  157°C 


17  hrs  108°C 
►24  hrs  158°C 


19  hrs  110°C 
►27  hrs  158°C 


19  hrs  110°C 
►27  hrs  158°C 


22  hrs  110°C 
►25  hrs  157°C 


22  hrs  110°C 
►25  hrs  157°C 


*  Cast  15  was  not  cast  under  vacuum  nor  was  vacuum  treatment  used 
after  casting.  The  remainder  were  either  poured  under  vacuum  or 
underwent  a  series  of  cycles  of  about  2  minutes  each  to  remove 


25 

0*058 

0*033 

0*028 

i  0*032 

0*030 

30 

0*082 

0*042 

0*042 

0*038 

0*036 

35 

0*125 

0*066 

0*064 

0*050 

0*048  j 

40 

0*220 

0*120 

0*100 

0*080 

0*088 

80 

0*122 

0*430 

0*190 

0*200 

85 

0*150 

0*102 

0*160 

0*090 

0*105 

1 

90 

0*090 

0*060 

.  0*070 

i  0*046 

0*060 

95 

0*052 

0*030 

0*036 

0*025 

0*030 

100 

0*033 

0*014 

0*021 

i 

0*013 

0*015 

_i _ _ _ 1 

,TP?  "T/tv’'®'!'  ”wrp r'ywT’-T!^ “ -' v’WiJ V’rW'yTfl, 7r"'*^r,» '  !T  "V, 


SPECIMEN 


CAST  No. 


TABLE  3 

MEASURED  VALUES  OF  FIBRE  CONCENTRATION 


METHOD 

USED 


Density 

Ash 

Mean 

Density 

Ash 

Mean 

Density 

Ash 

Mean 

Density 

Ash 

Mean 

Density 

Ash 

Mean 

Density 

Ash 

Mean 

Density 

Ash 

Mean 


Density 

Ash 

Mean 


Density 

Ash 

Mean 


Density 

Ash 

Mean 


Density 

Ash 

Mean 


SPECIMEN  ORIENTATION 
30°  1  45°  1  60c 

l  FIBRE  v/v 


3-15 

3*05 

3*05 

3-05 

3-02 

2-99 

2-93  j 

3-02 

3-1 

3*0 

3-0 

3*0 

2-65 

2*55 

2-55 

2-55 

2-58 

2-58 

2*59 

2-55 

2-6 

2-6 

2-6 

2*5 

8-7 

7*04 

8*70 

J  10*31 

8-34 

6-69 

8*52 

i  9*87 

8-5 

7-0 

8-6 

io-i 

12-68 

12-61 

12-7 

38*37 

38-4 

32-93 

32-80 

32-9 

41-27 

40- 95 

41- 1 


14-42 

14-29 

14-4 

38*67 

38-79 

38-7 

31-50 

31-24 

31*4 

36- 95 

37- 41 
37*2 


15*34 

15-02 

15-2 

41  43 
4U-  '6 
41-2 

31*50 

31-59 

31-5 

36*49 

36-5 


1  '  c2 
14  12 
14-5 

37-76 

37-8 

31-20 

31*41 

31-3 

37-26 

37*56 

37-4 


56- 43 
58-86 

57- 6 


53-87 

53-80 

53-8 


12-68 

12-67 

:,2-7 

t 

'4*59 
34-82  : 

34-7 

30- 28  ! 

31- 08 

30*7  | 

41-34  I 
44*98  i 


TABLE  3  (cont.) 


SPECIMEN 

%  FIBRE  v/v 

CAST  No. 

TEST  LENGTH 

METHOD  USED 

55  ran 

Density 

9-90 

17 

56  ran 

Density 

9-73 

40  ran 

Density 

11-63 

40  mm 

Density 

12-37 

40  mm 


TABLE  4 


EFFECT  Oi  SPECIMEN  LENGTH  ON  MEASURED 
AND  CALCULATED  SHEAR  MODULUS 


1  ■ 

» 

SPECIMEN 

TEMP. 

MEASURED 

CALCULATED 

CAST  No. 

TEST 

MODULUS 

MODULUS 

LENGTH 

c 

(GN/m2 ) 

(GN/m2  ) 

25 

1*31 

1-11 

55  mm 

100 

0-00421 

0*00298 

25 

1-27 

1-15 

25 

1-42 

1-11 

56  mm 

100 

0-00422 

0-00296 

17 

25 

1*39 

1-15 

25 

1-19 

1-15 

40  mm 

100 

O' 00415 

0-00308 

25 

0-96 

1-19 

25 

1-19 

1-17 

40  nm 

100 

O' 00407 

0-00313 

25 

0-99 

1-21 

25 

9-24 

- 

40  mm 

100 

0-00244 

- 

25 


35 

(Control) 


40  mm 


o  • 


TABLE  5 

CALCULATED  MODULUS  AND  TAN  6  VALUES  FOR  Su 
VARYING  FROM  J  S2 i  TO  2  Sj,  (CAST  20) 


TABLE  6 

EFFECT  OF  ASSUMING  DIFFERENT  VALUES  FOR 
POISSON'S  RATIO  ON  SHEAR  MODULUS  (G) 


POISSON'S 

RATIO 


TEST 

TEMP. 


SHEAR  MODULUS 
(GN/m1 ) 


ANGLE  BETWEEN  FIBRES  AND  SPECIMEN  AXES 
0°  I  30°  I  45°  |  60°  i  90° 


0»35 

0*40 

0*45 


45°C 


0*78 

0*78 

0*78 


1*08 

1*10 

1*12 


1-06 

1-08 

1*11 


0*86 

0*88 

0-91 


0-68 

0-69 

0*70 


0-35 

0-40 

0-45 


50°  C 


0*49 

0-49 

0-49 


0*7 

0*7 

0*7 


0-68 

0*69 


0-55 

0*56 

0-58 


0*43 

0-43 


I 

1 

3 


0*71 


0*44 


r 


TABLE  7 

TYPICAL  IMPACT  STRENGTHS  OF  VARIOUS  MATERIALS 


MATERIAL 


Stt.l 


IMPACT  STRENGTH 
(Nra/mm  WIDTH) 

IZOD  C HARPY 


Aluminium 

Phenol  formaldehyde  (unfilled) 
Phenol  formaldehyde  (fabric  filled) 
Phenol  formaldehyde  (laminates) 
Melamine 

Melamine  laminates 
Polyester 

Polyester  laminates 
Epoxy 

Epoxy  laminates 
Polyethylene  low  density 
Polyethylene  high  density 
Polyvinyl  chloride  rigid 
Polyvinyl  chloride  polyblend 
Nylon  66 

Nylon  (glass  filled) 

Polystyrene 
Polystyrene  (ABS) 

PMMA 


1-6 

to 

6-4 

1-1 

o-i 

to 

0-2 

0-016 

to 

0-02 

0-05 

to 

0-21 

0-54 

to 

r^» 

CM 

0*05 

to 

0-2 

0-02 

0*04  to  0*05 


0-01 

to 

0-02 

0-08 

to 

0*9 

0-01 

to 

0-05 

0-4 

to 

0-8 

0-86 

0-03 

to 

0-2 

0-02 

to 

0-05 

0-16 

to 

o 

CO 

O 

0-05 

0-16 

to 

CM 

CO 

o 

0-013  to  0-032 
0*05 

0-02  0-02 


1  J3 ^'ww^-'i'.'MvA^^J uw»j.u ; Jjiwi.f  ; IU  I  m'MM'UU.' JJ i 


TABLE  8 

TYPICAL  SURFACE  ENERGIES  OF  VARIOUS  MATERIALS 


MATERIAL 


Steel,  ship  plate  ductile 
Aluminium  alloy 

Polymethyl  methylacrylate  (tension) 
Polyesters 

Vulcanised  natural  rubber 
Glass  lantern  slide  covers  (moist) 
Cellulose  acetate 


CHARACTERISTIC  SURFACE 
ENERGY 

G0  (10"4.  J/W  ) 


3500 

260  -  1050 


8-7 

2-0 

120 

0*07 

50 


> 


k 


Polystyrene  cleavage 

Polymethyl  methacrylate  cleavage 
Polyester  cleavage 

Epoxy  cleavage 


4 

1*2 

0-1 


1 


Hot-stretched  polystyrene 
(parallel  to  stretch) 


Epoxy  (OER  332)  and  glass  laminate 
Impact  resistant  polystyrene  ABS 


Data  from: 


k  Irwin  (Ref. 35) 

1  McGarry  and  Self ridge  (Ref. 24) 
m  Broutman  (Ref. 25) 
n  McGarry  and  Self ridge  (Ref. 24) 


TABLE  9 


CALCULATED  VALUES  OF  6  [l-F(y)]J  f(y) 


n 

6[1-F(y)]*  f(y) 

—  00 

0 

-  5 

0-9  x  10"* 

-  4 

0*8  x  10"J 

-  3 

0*264  x  io_l 

-  2 

0*288 

-  1*5 

0*546 

-  1 

0*611 

-  0'5 

0*481 

0 

0*75  x  10“* 

1 

0*243  x  10-3 

FIG.  1  MOULD  FOR  UNIDIRECTIONAL  ARRAYS 


Reproduced  from 
best  available  copy. 


40 


FIG.  2 


FIG.  3 


FIG.  4 


41 


G,N/mz  TAN  4  G„w/m 


Hashin  —  Hill  bounds 


Fig  10.  Results  on  Specimens  cut  with  Axes  at  0°  to  Fibres 
Casts  13  to  20  Glassy  State 


Hashin-Hill  bounds 
Whitney— Riley  I 


Fig  11.  Results  on  Specimens  Cut  with  Axes  at  30  to  Fibres 
Casts  13  to  20  Glassy  State 


shin— Hill  bounds 
itney-Riley  I 


)  20  30  40  50  20  30  40  50  20  30  40  5  0  20  30  4  0  50  20  3  0  4  0  50 

CAST  13  3*/*  I  CAST  15  2-6*/.  I  CAST  16  1 0*1  •/•  I  CAST  18  12*7»/.  CAST  20  34*7 

Temperature  *C 

Fig  14.  Results  on  Specimens  Cut  with  Axes  at  90°  to  the  Fibres 

Casts  13  to  20  Glassy  State 


17.  Results  on  Specimens  Cut  with  Axes  at  45°  to  Fibres  modutul  cqlrumt 

Compared  with  Theoretical  Calculations  and  Control  from  thret  tfteor 

Casts  13  to  20  Rubbery  State 


Fig  IB.  Rfiulh  on  SpfcIrMM  Cur  »llh  Ahs  nl  du  f»  Fibres  fc»l\W  fflOdwIul  cptculoled 

C«mp?r«d  wtHi  T^twnricflf  Cfkutetim  and  Confi^l  from  Iftrw  theone* 

Cqh  13  lo  20  Rubbaiy  Stott 


10  20  30  40  50  60  70  80 

Temperature  *C 


90  100  110  120 


ig  20.  Comparison  of  results  on  specimens  cut  with  Axes  at  0°  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Glass 
F  ibres 


Cast  29 ,  33%  F  ibres 


10  20  30  40  30  60  70  80  9  0  100  110  120 

Temperature  *C 


Fig  21  Comparison  of  results  on  Specimens  cut  with  Axes  at  30°  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Glass 
Fibres 


Cast  29  31%  Fibres 


10  20  30  40  50  60  70  80  90  100  110  120 

Temperature 


Fig  22  Comparison  of  results  on  Specimens  cut  with  Axes  at  45°  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Glass 
Fibres 


Cast  29  32%  F  ibres 


10  20  30  40  50  60  70  80  90  100  110  120 

Temperature  *C 


Fig  23  Comparison  of  results  on  Specimens  cut  with  Axes  at  60°  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Glass 
Fibres 

Cast  29  31%  Fibres 


10  20  30  40  56  60  70  80  56  100  110 

Temperature  *C 


Fig  24.  Comparison  of  results  on  Specimens  cut  with  Axes  at  90°  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Fibres 


20  30  40  90  60  70  90  90  100  110 

Temperature  *C 


25.  Comparison  of  results  on  Specimens  cut  with  Axes  at  0°  to  Fibres 
with  Theoretical  Calculations  and  Control-Finish  on  Fibres 


X)  20  30  4090  to  70  80  90  100  110  ItO 

Temptfaturt  *C 


Fig  26.  Comparison  of  rosults  on  Specimens  cut  with  Axes  at  45°  to  Fibres 
with  Theoretical  Calculations  and  Control -Finish  on  Fibres 
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Fig  27,  Comparison  of  results  on  Specimens  cut  with  Axes  at  0°  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Fibres 
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Fig  28 


Comparison  of  results  on  Specimens  cut  with  Axes  at  0  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Fibres 


Cast  42  58%  Fibres 
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Comparison  of  results  on  Specimens  cut  with  Axes  at  0°  to  Fibres 
with  Theoretical  Calculat  ions  and  Control  -  Finish  on  Fibres 


Cost  34  43‘X  Fibres 


o  0  °  \  ° 

\ 

\ 

Control  -  n 

Experiment  o,x - 

heory  n 

-> - 1 - »--  ■  j  i  i  i - j, - 1— ...  i. _ i 

20  30  40  30  60  70  80  90  100  110  120 

Temperature  *C 


Comparison  of  results  on  Specimens  cut  with  Axes  at  0°  to  Fibres 
with  Theoretical  Calculations  and  Control  -  No  Finish  on  Fibres 


rf  =  lOOji 
£  =  0.01 


Fig  32.  Graph  of  Predicted  Failing  Strain  against  Fibre  Concentration  for  rf  =  lOp-,  £=0.01 


Fig  34.  Graph  of  Predicted  Failing  Strain  against  Fibre  concentration  for  rf  -  10p.» 


